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We present a theoretical model which accounts for recent angle resolved photoemission data in 
bilayer cuprate superconductors. Lineshapes and dispersions of the various bonding and antibond- 
ing features in the spectra are quantitatively reproduced. The observed dispersion anomalies are 
consistent with the interaction of electrons with a bosonic mode which is odd with respect to the 
layer indices, a unique property of the magnetic resonance observed by inelastic neutron scattering. 
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Recent angle resolved photoemission (ARPES) exper- 
iments on bilayer cuprate superconductors were able for 
the first time to resolve a bilayer splitting between bond- 
ing and antibonding bands ||[ |[ [|. The dispersion 
near the (n, 0) point of the Brillouin zone shows an un- 
usual asymmetry between bonding and antibonding self 
energy effects. In particular, Feng et al. [Q found that 
the energy distribution curves (EDCs, taken at constant 
momentum) in overdoped Bi 2 Sr 2 CaCu208+5 {T c — 65K) 
consist of three features: an antibonding band (AB) 
peak near 20 meV, a bonding band (BB) peak near 
40 meV, and a bonding hump near 105 meV. Recently, 
Gromko et al. Q reported strong self energy effects in the 
dispersions derived from momentum distribution curves 
(MDCs, taken at constant energy) in similar samples 
(T c = 58K). Near momentum (k x ,k y ) = (1, 0.13)7r/a, 
an ^-shaped dispersion anomaly, discussed previously in 
Ref. Q, was shown to be present only in the bonding 
band MDC, at binding energies between 40 meV and 
60 meV. In both experiments, a low energy double peak 
structure in the EDC was only resolvable in the same 
momentum region. 

In this letter, we demonstrate that all of these fea- 
tures can be explained by a model which assumes that 
low energy scattering of electrons between the bonding 
and antibonding bands is strong compared to scattering 
within each of those bands. As scattering events which 
connect different bilayer bands are odd with respect to 
permutation of the layers within a bilayer, this implies 
that the corresponding bosonic excitations which medi- 
ate such scattering must be dominant in the odd chan- 
nel. The magnetic resonance observed in inelastic neu- 
tron scattering has exactly this property || [tJ] , and more- 
over has the correct energy to reproduce the observed 
dispersion anomalies. Thus, the recent ARPES experi- 
ments on bilayer cuprates provide independent support 
for a strong coupling between this resonance and elec- 
tronic excitations. 

If the electrons are phase coherent between the two 
planes, then the spectra will exhibit separate bonding (b) 
and antibonding (a) features with (normal state) disper- 
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In the superconducting state, the dispersions are modi- 
fied by the presence of the d-wave order parameter given 
by Ag = Ao(cosfca; — cosfcy). In agreement with ex- 
periment [3|, we assume Ao to be the same for the 
bonding and antibonding bands. Then, the dispersion 
in the superconducting state takes the form (ijl ' 6 -*) 2 = 

i&y + (A_g) 2 . 

This "non-interacting" picture, though, is insufficient 
to describe the three observed dispersion features. We 
are able to account for them in a model where electrons 
are coupled to a resonant spin mode situated below a 
gapped spin fluctuation continuum. Such a spectrum is 
observed in inelastic neutron scattering experiments ||. 
In bilayer materials, the spin susceptibility is a matrix in 
the layer indices, having elements diagonal {xaa, Xbb) and 
off-diagonal (x&a, Xab) in the bonding- antibonding rep- 
resentation. The components of the spin susceptibility 
transforming as even and odd with respect to the layer 
indices are given by Xe = Xaa + Xbb and Xo = Xab + Xba- 
For identical planes, Xaa = Xbb and Xab = Xba- The 
measured susceptibility is then given by 



X = Xe COS 
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where d is the separation of the layers within a bilayer. 
The resonance part, Xres, was found to exist only in the 
odd channel, whereas the continuum part, Xc enters in 
both ||. Thus, Xo = Xres + Xc and Xe = Xc- This 
means the resonance mode can only scatter electrons be- 
tween the bonding and antibonding bands. In contrast, 
the spin fluctuation continuum scatters both within and 
between these bands. As we demonstrate below, the odd 
symmetry of the resonance is crucial in reproducing the 
ARPES spectra. 

We employ for the single layer self energy the func- 
tional form of Ref. B . It is proportional to the convolu- 
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tion of the Gor'kov-Green's function G with the dynamic 
spin susceptibility x- We write this self energy symboli- 
cally as S = g 2 x*G (the hat denotes the 2x2 particle hole 
space, and g is the coupling constant). For a bilayer sys- 
tem, there is a separate self energy for each band, & a ' b \ 
Using the above notation, these self energies are given by 

= 9 Y{xre S *G^+ X c*(G^+G^)} (3a) 

± (a) = ^{x,e S *GW+Xc*(G (6) +G (a) )}.(3b) 

Dispersion anomalies arise mainly from coupling to the 
resonance mode. This means that dispersion anomalies 
in the bonding band are determined by the antibond- 
ing spectral function and vice versa. (The spectral func- 
tions for the bonding and antibonding bands are given 
by the imaginary parts of the diagonal components of 
the Green's functions.) Because the antibonding band is 
(in contrast to the bonding band) close to the chemical 
potential at (n, 0) [0], the associated van Hove singularity 
leads to a larger self energy for the bonding band. 

The self energy effects can be described very accu- 
rately by using a model spin spectrum which consists 
of a perfectly sharp resonance mode below a gapped con- 
tinuum with gap Ah- We approximate the intensity of 
the continuum to be constant. In reality, the intensity 
of the continuum decays at high energies. However, as 
we are interested in self energy effects in an energy range 
|e| < 200 meV, the above approximation is adequate. A 
high frequency cut-off in the convolution integrals in Eqs. 
( |3a]) and (|3t]) was introduced. The precise cut-off proce- 
dure however does not affect the low energy physics. Any 
variation in the cut-off can be accounted for by a read- 
justment of the coupling constants and the normal state 
renormalization factor. We model the resonance mode, 
R = 2Imx reS ) and the continuum, C — 2Imx c , by 

Rlu,q = 2r(q){S(u -Sl res ) - 5(-oj - Q res )} , (4a) 
C W|? = 2c(q) {6(w - 2A h ) - G{-w - 2A h )} . (4b) 

where the function is zero for negative argument and 
one otherwise. The momentum dependences of the res- 
onance mode and the continuum, given by the functions 
r(q) and c(q), are plotted in Fig. 0. The resonance mode, 
shown in Fig. 0a, is peaked at Q, with a correlation length 
of £res — 2a, where a is the lattice constant. The gapped 
continuum, shown in Fig. 0b, is much broader. This is 
motivated by the experimental data |3| , which show that 
the continuum is enhanced around Q with a correlation 
length of only 0.5 lattice constants. Also, the momentum 
dependence of the continuum excitations exhibit experi- 
mentally a flat behavior around Q, as in Fig. |b. To sim- 
plify the model, we use the same functional form for the 
gapped continuum in the even and odd scattering chan- 
nels. This is consistent with the superconducting state 
data, where the gap in the odd channel (about twice the 




FIG. 1: Momentum dependence of a) the resonance mode 
and b) the gapped spin fluctuation continuum. The resonance 
mode is peaked at Q = (n, ir) with a correlation length equal 
to twice the lattice constant, £ res = 2a. The continuum spec- 
trum, in contrast, is rather broad around Q. 

maximal superconducting gap), is close to the optical gap 
in the even channel. 

We have chosen parameters appropriate for the over- 
doped sample (T c — 65K) studied in experiment Q. The 
normal state dispersion is obtained from a 6 parameter 
tight-binding fit to experimental data jflj, with the bi- 
layer splitting given by Eq. El In order to determine the 
total of 7 parameters for this fit, we used three Fermi 
wavevectors, k N = (0.37, 0.37)7r/a, k^ } = (1, 0.217)7r/a, 
and = (1, 0.135)7r/a, the Fermi velocity at k^, 

vn = 2.316 eVA, and the normal state band energies 
at (tt,0), Cm = -105 meV and = -18 meV (the 
corresponding bilayer splitting is — = t± = 87 
meV). Finally, to obtain the correct overall shape of the 
dispersion, we fixed the band energy at (it, 7r), using a 
reasonable value £y = 0.8 eV. The high energy (|e| > 200 
meV) dispersions are not affected strongly when going 
from the normal to the superconducting state. How- 
ever, even in the normal state, the bare dispersion is 
renormalized by the normal state spin fluctuation contin- 
uum. To account for this renormalization, we multiply 
the above dispersion by a factor of 1.4 in order to ob- 
tain the bare dispersion. For the remaining parameters 
of the model, we use Ao = 16 meV Q, £! res = 27 meV, 
A h = 0.9A , g 2 r(Q) = 0.15eV 2 , and g 2 c(Q) = 0.72eV. 
The value for the resonance energy was obtained from the 
relation Q res — 4.9fcsT c found experimentally to hold for 
overdoped Bi 2 Sr2CaCu208+5 §, 0, [T^]. The resonance 
weight r(Q) has not been measured for overdoped materi- 
als. For optimally doped Bi 2 Sr 2 CaCu208+(5, it is 0.95^ 
per plane 0, and we estimated for this case g = 0.65 
cV (9|. Assuming for the overdoped sample the same 
coupling constant, the above value for g 2 r(Q) would im- 
ply a resonance weight of 0.36/ig per plane. Similarly, 
with this coupling constant and our value for g 2 c(Q), we 
obtain a (2D)-momentum averaged continuum contribu- 
tion of 1.7/j,^/eV per plane (gotten by summing the even 
and odd channels for energies u> <0.2eV). Our calcula- 
tions are for a temperature T = 10K. Note we use un- 
renormalized Green's functions in Eqs. Ba] and Bn. This 
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FIG. 2: a) Tight binding Fermi surfaces for antibonding and 
bonding bands in overdoped Bi2Sr2CaCu2C>8+<5 (T c = 65K). 
b) Renormalization function and c) imaginary part of the self 
energy at the (n, 0) point of the zone for the bonding band 
(BB) and antibonding band (AB). d) Spectral functions for 
k y — n/a, and k x varying from to 0.277r/a. For comparison 
with experiment, the spectra are convolved with a Lorentzian 
energy resolution function (FWHM 12 meV). 
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FIG. 3: Dispersion of EDC peak positions (left) and MDC 
peak positions (right) near the (n, 0) point of the Brillouin 
zone. The EDC dispersion consists of three branches, one 
antibonding peak, one bonding peak and one bonding hump. 
The MDC bonding dispersion shows a characteristic S shape 
behavior. The small symbols in the right picture show again 
the EDC peak positions for convenience. 



approximation is sufficient to explain a large variety of 
data, and can be justified by considerations discussed in 
Ref. 0. 

The bonding and antibonding normal state Fermi sur- 
faces are shown in Fig. ||a. The bilayer splitting is max- 
imal near the (tt, 0) points of the zone. Thus, we will 
concentrate on this region in the following. In Fig. |^b we 
show the renormalization function for bonding and anti- 
bonding bands at the (tt, 0) point. The renormalization 
is stronger for the bonding band than the antibonding 
band. This is a result of the proximity of the antibond- 
ing saddle point singularity to the chemical potential. As 
is seen in this figure as well, both bands are renormalizcd 
up to high energies. The high energy renormalization ap- 
proaches that of the normal state dispersion (1.4). 

The imaginary part of the self energy is shown in 
Fig. ^c for the bonding and antibonding bands. As emis- 
sion processes are forbidden for |e| < fl re s, the imaginary 
part of the self energy is zero in this range. Due to scat- 
tering events to the antibonding band, electrons in the 
bonding band have a large imaginary part of the self en- 
ergy in the range between 40 and 60 meV. These events 
are dominated by emission of the resonance, and are en- 
hanced due to the van Hove singularity in the antibond- 
ing band close to the chemical potential. In contrast, the 
imaginary part of the antibonding self energy is not en- 
hanced because the bonding band is far away from the 
chemical potential at (tt,0). Consequently, it shows lin- 
ear behavior over a wide energy range, with a gap at low 



energies (|e| < fl re s)- 

Fig. ||d presents the intensities for the bonding and 
antibonding spectra. The antibonding spectra consist of 
a low energy AB peak, and the bonding spectra have 
a low energy BB peak and a higher energy BB hump 
feature. In agreement with experiment ([fij and j|]), the 
width of the EDC spectrum is large for the BB hump, but 
not so for the BB and AB peaks. We also mention that 
the BB peak is well defined in the whole region between 
the BB Fermi crossings on either side of (tt, 0), but the 
finite energy resolution does not allow to resolve it near 
(tt, 0) as seen in Fig. ||d. 

In Fig. H we show our results for dispersion of the 
EDC and MDC peak positions. In the EDC disper- 
sions, Fig. |j|a, we reproduce the experimentally observed 
three branches jl| , one antibonding peak and two bond- 
ing branches, denoted 'BB peak' and 'BB hump'. The BB 
peak has a very flat dispersion near k x = in agreement 
with experiment Q. Its position at 40 meV is approxi- 
mately given by Sl r es + A^, where is the gap at the 
antibonding Fermi crossing. Thus, the energy separation 
between the AB peak at the AB Fermi crossing and the 
BB peak at (tt, 0) is a measure of the resonance mode en- 
ergy fl res m overdoped compounds. The BB hump posi- 
tion at high binding energies (105meV) is determined by 
the normal state dispersion of the bonding band. Because 
the spin fluctuation continuum changes only at low ener- 
gies when going from the normal to the superconducting 
state, the position of the BB hump maximum is not very 
different from the normal state BB dispersion. This is in 
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FIG. 4: Spectral functions for k y — Tv/a at k x — 0, 
fca: = 0.127r/a and fen = 0.2l7r/a. The full lines are the sum 
of antibonding (dotted) and bonding (dashed) spectral func- 
tions. For comparison with experiment, a Lorentzian energy 
resolution function of 12 meV width was assumed. The dou- 
ble peak structure is clearly resolved for k x = 0.12n/a, as in 
experiment ^, B|. 



agreement with experiment . The intensity of the AB 
peak decreases quickly when it approaches the BB peak, 
but is strong at (tt, 0) because of the proximity of the AB 
band to the chemical potential in this region. 

In Fig. ||b, we present the MDC dispersions (for com- 
parison we also reproduce the EDC dispersions as small 
symbols). The MDC dispersion consists of two branches, 
an AB MDC branch and a BB MDC branch. The self on 



indicated. The spectra are convolved with a Lorentzian 
energy resolution function to allow for direct comparison 
with experiment. We reproduce all experimental find- 
ings. First, at (ir,0), only the BB hump and the AB 
peak are resolved. This is due to resolution effects men- 
tioned above. Second, near the AB Fermi crossing, the 
spectra show a characteristic double peak structure, with 
a relatively sharp AB peak and a BB peak separated from 
a broad BB hump. Third, at the BB Fermi crossing, only 
the BB peak is observed. The BB hump is so small in 
intensity that it only leads to a kink-like feature in the 
spectrum. 

We have presented a theory to account for the exper- 
imentally observed self energy effects in the bilayer split 
bands in double layer high temperature superconductors. 
We reproduced quantitatively the EDC dispersions, the 
MDC dispersions, and the spectral lineshapes. We found 
that the ARPES data are consistent with the interaction 
of the electrons with a sharp bosonic mode which is odd 
in the layer indices, a property unique to the magnetic 
resonance observed by inelastic neutron scattering. 
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